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Abstract. We consider Voronoi's reduction Aeory of positive definite quadratic 
forms which is based on Delone subdivision. We extend it to forms and Delone 
subdivisions having a prescribed symmetry group. Even more general, the theory is 
developed for forms which are restricted to a linear subspace in the space of quadratic 
forms. We apply the new theory to complete the classification of totally real thin 
algebraic number fields which was recently initiated by Bayer-Fluckiger and Nebe. 
Moreover, we apply it to construct new best known sphere coverings in dimensions 
9,. ..,15. 
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1. Introduction 

In this paper we generalize a classical reduction theory for positive definite qua- 
dratic forms due to Voronoi BVorOSi His theory gives in particular an algorithm to 
classify Delone subdivisions of R"^ with vertex-set Z*^ up to the action of GLrf(Z). For 
precise definitions of used terms and a brief description of the classical theory we refer 
to Section 121 

We present our generalization in Section |3] We extend the classical theory in two 
different directions. On the one hand, we generalize the theory from vertex-set to 
general periodic vertex-sets. On the other hand we give an equivariant theory dealing 
with positive definite quadratic forms with a prescribed automorphism group G < 
GLf;(Z). This equivariant theory is an analogue of the theory of G-perfect forms by 
A.M. Berge, J. Martinet and F. Sigrist, [.BMS92,1 . which was motivated by the search of 
good packing lattices with prescribed symmetries. In fact, as in their case, our theory 
can be developed in the more general context of a linear subspace T of quadratic forms. 

As in the classical theory, where T is the space of all quadratic forms, we get a 
polyhedral subdivision of the space of positive definite forms in T into generic T- 
secondary cones (in the classical case also called L-type domains), which contain those 
forms giving the same Delone subdivision. In contrast to the classical theory these 
generic subdivisions ai^e no longer triangulations. In the equivariant theory, we have 
only finitely many T-secondary cones up to the action of GLrf(Z). Our proof of this 
fact in Section |4] uses only the action of GLd(Z) on a polyhedral subdivision of the 
space of positive definite quadratic forms. So it applies to the theory of G-perfect 
forms and it gives a unified view on both theories (see Remark l4!4l ). 

We describe the theory in a way which allows us to work with it computation- 
ally. In particular we made some effort to reduce redundancies in the description of 
secondary cones (see Theorem 13.11 ). Also the transition from a Delone subdivision 
of a T-secondaiy cone to the Delone subdivision of a contiguous T-secondary cone, 
called a T-flip, is given explicitly (see Theorem 13.81 ). For it we define repaitition- 
ing polytopes, in which a polyhedral subdivision has to be replaced by another. As 
a nice byproduct we obtain an explicit description of fiips which occur in the theory 
of equivariant secondary polytopes of regular subdivisions (of a polytope), recently 
introduced by Reiner IIRei02l (see Remark [l!9l ). 

We are not the first who consider generalizations of Voronoi 's theory. Periodic 
tilings and reduction theory of positive definite quadratic forms appear naturally in 
algebraic geometry in the study of degenerations of abelian varieties and compact- 
ifications of Siegel modular varieties. For these reasons, Delone subdivisions and 
Voronoi 's reduction theory were used and studied by Mumford and Namikawa about 
25 years ago, and by many algebraic geometers since then. In Section [5] we review 
generalizations which came up in this context and compare them to our work. 

We apply our extension of Voronoi's theory to two different problems. We use it to 
finish the classification of totally real thin number fields, which was recently started 
by Nebe and Bayer-Fluckiger IIBN05II (Section [H). We use it to construct new best 
known sphere coverings (Section |9l). Both applications involve finding best lattice 
configurations with respect to a given Delone subdivision. A brief description of the 
problem and of the convex optimization tools we used is given in Section |6] Since 
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we are dealing with non-linear optimization problems we can usually only approxi- 
mate the lattices we want to find. By using convex optimization duality and rational 
approximations we can give mathematical rigorous error bounds for the quality of ap- 
proximated lattices. Algorithmic issues in the classification of T-secondary cones are 
addressed in Section |7i 

So far, using the new theory for the lattice case, we found new best known lattice 
sphere coverings in dimensions d = 9, . . . , 15. Using the classical theory (cf. IISV06II ) 
and new methods to enumerate all vertices of symmetric Voronoi cells efficiently (cf. 
IIDSV07II ). we found new best known coverings in dimensions d = 6, 7, 8 and d = 
17, 19, 20, 21 as well. A complete list of the best known values is given in Table 2 in 
Section |9] With the exception of dimension d = Q and d = 7, we do not think that 
these lattice coverings are optimal. 

Furthermore, we strongly believe that our extension of Voronoi 's classical theory to 
periodic sets is a first step towards the construction of non-lattice coverings, which are 
less dense than any lattice covering. 

2. Background: Lattices, PQFs and Delone subdivisions 

We start with basic definitions and basic results and some background on Voronoi's 
reduction theory. In the first section we introduce lattices and positive definite qua- 
dratic forms, PQFs from now on. For further reading we refer to IICS88II and IIGL87II . 
In the second section we introduce Delone polyhedra and Delone subdivisions. For 
further reading about concepts related to polyhedra we refer to IIZie97i In the third 
section we give a very brief review on Voronoi's results in irVbrOSi 

Let be the d-dimensional Euclidean space with column vectors x = (xi , . . . , x^)* 
and norm ||a;|| = Vxhc. 

2.1. Lattices and PQFs. A d-dimensional lattice L in is a discrete subgroup L = 
TLv\ + • • • + ZtJrf with linearly independent Vi G W^. The family (ui, . . . is 
called a basis of L. To it we associate the positive definite symmetric Gram matrix 
Qb = B^B, where B G GL£;(M) is the invertible matrix whose i-th column is Vi, and 
L = BZ'^. 

Given on the other hand a positive definite symmetric matrix Q, there exists a matrix 
B G GLci(M) with Q = B^B. The matrix B is uniquely determined up to orthogonal 
transformations. Any other A G GLrf(IR) with AZ'^ = BU^ can be written a.?, A = BU 
with U G GLrf(Z). This relation yields Qa = A* A = U^QbU. We say Qa and Qb 
are arithmetical equivalent in this case. 

The space of real symmetric d x d matrices is denoted by S^. It is a C^^^)- 
dimensional Euclidean space with inner product {A,B) = trace(Ai?). The subset 
SyQ of positive definite symmetric matrices is an open convex cone in S"^. Abusing 
notation, we identify positive definite quadratic forms, PQFs from now on, and positive 
definite symmetric matrices by 

Q[x] = x^Qx = {Q,xx^). 

The topological closure of S^q in S'^ is the set 5>q of all positive semidefinite matri- 
ces. By the relations above, 5>q can be identified with Od(]R)\GLfi(]R), where Orf(M) 
denotes the subgroup of orthogonal dx d matrices in GL(i(IR). The group GLrf(Z) acts 
on SyQ by Q 1-^ U^QU. Thus the set of isometry classes of d-dimensional lattices 
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(i.e. (i-dimensional lattices up to orthogonal transformations) can be identified with 

2.2. Delone polyhedra and Delone subdivisions. A polyhedron is a set in which 
can be represented as a finite intersection of closed half spaces, e.g. given by a system 
of linear inequalities. A polytope is the convex hull of finitely many points and by 
a theorem due to Minkowski and Weyl, polytopes are bounded polyhedra and vice 
versa. Given a discrete set A C W^, a polytope P = convji?!, . . . , Vn} with vertices 
t;j G A is called a Delone polytope, if there exists a center c and a radius r > 
such that the Euclidean distance between c to all points ?; G A satisfies ||c — ?;|| > r, 
with equality only for the vertices of P. The set of all Delone polytopes of a lattice 
forms a polyhedral subdivision of R'^. This is a family of polyhedra, called faces, 
whose union is and which is closed with respect to intersections. Each face of the 
subdivision sharing relative interior points with another face of its dimension coincides 
with this face. Note, by a theorem of Gruber and Ryshkov IIGR89II . the latter face-to- 
face property holds whenever it holds for the facets (faces of co-dimension 1). 

Our main interest is in vertex-sets A^ which are periodic, that is, a finite union of 
lattice translates of a lattice L, e.g., A^, = IJI^i U + L with tj G R"' for i = 1, . . . , m. 
In many cases it is convenient to work with coordinates with respect to a given basis 
of L, say given hy B ^ GL(i(M). This means that we work with standard periodic sets 

m 

(1) A = y + 

i=l 

and norm defined hy Qb- Hence the norm of a; G M'' is given by \/Qb[x]- A polytope 
P = conv{Di, . . . , Vn} with Vi G A, is called uDelone polytope of Qb and A, if there 
exists a c G M'^ and a real number r with Qb[c — v] > r^, where equality holds if and 
only if 17 is a vertex of P. 

It is a bit more general (and in some situations convenient) to consider Delone sub- 
divisions of positive semidefinite forms Q as proposed by Namikawa in IINam76ll . Pos- 
itive semidefinite forms define seminorms on M'^ by a/ Q[x]. A (possibly unbounded) 
polyhedron P = conv{t7i, 1^2, • • •} with Dj G A is then called a Delone polyhedron of 
Q, if there exists a c G M'^ and a real number r with Q[c — Vi] = for i = 1, 2, . . . 
and Q[c — v] > r"^ for all d G A \ {vi,V2, ■ ■ ■}■ Suppose for Q G 5>q there exists a 
matrix U G GLrf(Z) and aQ' £ S^q with 

(2) U^QU = (0 = Q. 

Then d' = rankQ and if P is a Delone polyhedron of Q then U^^P is a Delone 
polyhedron of Q. Latter are of the form M''"'^ x P' with P' £ a Delone polytope 
of Q'. The set of all forms being arithmetical equivalent to some positive semidefinite 
form Q of the form Q, with Q' positive definite, is called the rational closure of 
SyQ, denoted by S^q. The importance of the rational closure is due to the following 
proposition. 

Proposition 2.1. For Q G 5>q exists a d-dimensional Delone polyhedron with respect 
to a standard periodic vertex-set if and only if Q £ 5>q. 
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Although this proposition might be known we are not aware of a reference. In 
IINam76[ §2.1] Namikawa showed that every Q lying in the rational closure has a d- 
dimensional Delone polyhedron. For completeness we repeat his argument. 

Proof. If Q £ 5>Q, there exist d-dimensional Delone polyhedra which after a suitable 
transformation in Q>Ld{1) are those of the form Q in Q. If on the other hand Q 5>q, 
then for all arithmetical equivalent forms Q as in (O, the form Q' £ 5>'q is not positive 
definite, hence d! > rank((5). For an aiithmetical equivalent form Q with minimal d' 
there exists no rational, hence no integral vector in the kernel of Q' . As a consequence, 
for such a Q' we find for all c G M'^ and all r > a i> G with Q'[v — c\ < r (see 
IISie89[ Lecture VI]). Therefore there do not exist -dimensional Delone polyhedra 
for Q' and Z'^ , and hence no d-dimensional ones for Q and Z'^ as well. The same is 
true for Q and a standard periodic set A. □ 

The set 5>q can also be described as the set of all non-negative combinations (the 
cone) of rank-1 forms vv^ with v ^ Z'^: 

Proposition 2.2. We have 



Proof. By the definition of 5>q, every Q G 5>q is arithmetical equivalent to a form Q 
as in dJll. The PQF Q' is a sum of rank-1 forms (cf. IIVor07[ Section 24]). Therefore Q 
and Q are of this form as well. 

Suppose on the other hand that Q G 5>q is the sum of rank 1-forms, e.g. Q = 
Yl^i o^i'^i'^i with Vi G and Oj > for z = 1, . . . , m. If rankQ = d there is 
nothing to show. If rank Q < d, then there exist linearly independent Pi, . . . ,pi, € 
with k = d — rankQ such that Q[pj] = Yl^i '^iip'jVi)'^ = for j = 1, . . . ,k. 
Thus the pj are orthogonal to each of the Vi and therefore they span a /c-dimensional 
linear subspace which contains a /c-dimensional sublattice of Z'^. We choose a basis 
(til, ■ ■ ■ , Uk) of this sublattice and extend it to a basis U = (iti, . . . , Ud) of U^. Then 
(C/-i)*Q([/-i) is of the desired form ©. □ 

The set Del((5) of all Delone polyhedra of a Q G 5>q is called the Delone subdi- 
vision of Q. If all elements in Del((5) are simplices, DeT((5) is called a triangulation. 
The subdivision Del((5) is a polyhedral subdivision of which is invariant under 
translations of the form x ^ x + v, where v G ifi. Therefore T>e\{Q) is completely 
determined by the stars of the translation vertices in ([U for i = 1, . . . , m, where 
a star of a single vertex is the set of all Delone polyhedra containing it. We call two 
Delone polyhedra P and P' equivalent if there exists a i> G so that P = v + P' . We 
say that Del((5) is a refinement of De\{Q') (and Y)e\{Q') is a coarsening of T)e\{Q)), if 
every Delone polytope of Q' is contained in a Delone polytope of Q. In Section [7] we 
need an algorithm which computes the Delone subdivision of a given PQF. The inter- 
ested reader can find a discussion of these computational issues in our paper IIDSV07II . 

2.3. Voronoi's reduction theory. Before we generalize Voronoi's reduction theory 
in the next section, we briefly recall the original theory (see HVorOSII . IIDel37ll and 
IISV06II ). Generally, the task of reduction is to find a fundamental domain in with 
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respect to the action of GLd{Z). Voronoi's reduction is based on secondary cones, also 
called L-type domains, of Delone triangulations with vertex-set Z*^. More generally, 
the secondary cone A('D) of a Delone subdivision V with a standard periodic vertex- 
set is defined by 

A{V) ={Qe Sio ■■ Del(Q) = V}. 

We say that two secondary cones of Delone triangulations are bistellar neighbors if 
the Delone triangulations differ by a bistellar flip, which is a specific change of the 
triangulation (see Section [33] for a definition and generalization). Voronoi also showed 
that the topological closures A('D) of secondary cones of Delone triangulations form 
a polyhedral subdivision of 5>q. 

Theorem 2.3 {Voronoi's Reduction Theory). 

The secondary cone of a Delone triangulation with vertex-set Z'^ is a full-dimensional, 
open polyhedral cone in S^q. The topological closures A(P) give a polyhedral subdi- 
vision ofS^Q. The closures of two secondary cones have a common facet if and only if 

they are bistellar neighbors. The group GLrf(Z) acts on the tiling by U A{T))U. 
Under this group action there are only finitely many inequivalent secondary cones. 

Note that by Voronoi's theory we have a non-intersecting subdivision of 5^o> as 
well as of 5>Q, into secondary cones. In it, every cone is an open polyhedral cone 
with respect to its affine hull. We refer to such a decomposition of S^q as an open 
polyhedral subdivision. Such subdivisions are of particular interest in Section |4l if 
they fall into only finitely many orbits under the action of GLrf(Z), as in the case of 
secondary cones. 

3. Generalization of Voronoi's reduction theory 

In this section we generalize Voronoi's reduction theory. For our generalization 
we consider a linear subspace T C S'^ and look at T-secondary cones of Delone 
subdivisions V defined by 

At{V) = A(P) n T. 

We call a T-secondary cone Ar(D) and the corresponding Delone subdivision V T- 
generic if dim Ay ("D) = dimT. By Voronoi's Theorem 12. 31 the topological closures 
of T-generic, T-secondary cones give a polyhedral subdivision of 5>q H T. Two T- 
generic cones ai^e called contiguous if their closures share a facet. difference with 
the classical theory is the existence of dead-ends, which are facets only incident to one 
T-generic cone. These necessarily contain only non-positive forms in 5>q \ S^q. 

The ultimate goal would be to state for every subspace T a theorem as Theorem 
12.31 which deals with the case T = S'^. It turns out though that in general, this is not 
always possible. If dim T = 1, the intersection of T with S^q contains a PQF Q and 
all its multiples. In this case a generalized Theorem 12. 3 1 is trivially true. Therefore, if 
not stated otherwise we assume dimT > 2 in what follows. 

The "road map" for our generalization is the following: In Section [3TT] and ll!2] we 
determine the secondary cone of an arbitrary Delone subdivision explicitly, because in 
our more general setup we have to deal with Delone subdivisions which are not Delone 
triangulations. Then in Section [331 we generalize the notion of bistellai^ neighbors to 
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the new setting. In the new theory not every subspace T gives a polyhedral subdivi- 
sion of 5>Q n T with only finitely many inequivalent T-secondary cones (see Remark 
14. lb . For specific vertex-sets and subspaces though, there exist only finitely many in- 
equivalent T-secondary cones. Such a finiteness result is given in Section|4]for and 
subspaces containing all PQFs which are invariant under a given finite subgroup of 
GLrf(Z). Thus we obtain an equivariant version of Voronoi's reduction theory. Using 
modern terminology, our proofs not only generalize Voronoi's theory, but also shorten 
his argumentation. 

3.1. Polyhedral description of secondary cones. Let P be a Delone subdivision 
with a standard periodic vertex-set A. In this section we want to describe the sec- 
ondary cone A('D). It will turn out that A (2?) forms a (relative) open polyhedral 
cone in Si^. Theorem O gives the precise statement. This result is an adaption 
of Voronoi's "fundamental theorem" (see HVorOSI §77]), which deals with the generic 
case of Delone triangulations with respect to the vertex-set Z'^ (see also IISV061 Section 
5.1]). Actually the first statement of the second part of Theorem 13. II is not explicitly 
stated in Voronoi's work. It goes back to Nakamura IINak75[ Lemma 1.1]. 

We describe below the polyhedral cone explicitly by Unear equalities and inequal- 
ities. The description is needed for our application and therefore we put some effort 
into avoiding redundancies. The linear equalities are coming from d-dimensional non- 
simplicial polyhedra in V. Hence, in the generic case of Delone triangulations there 
are no linear equalities. The linear inequalities are coming from {d — 1) -dimensional 
polyhedra (facets) in V. For the formulation of these linear conditions we define for an 
affinely independent set 1/ C M*^ of cardinality d+\ and a point w the quadratic 
form 

(3) Nv,w = WW* - ^ avvv\ 

where the coefficients a„ are uniquely determined by the affine dependency w = 

T,vev "-""^ with 1 = E^GV 

The following theorem generalizes Voronoi's "fundamental theorem" for Delone 
triangulations to arbitrary polyhedral subdivisions. As vertex-sets A C M'^ we allow 
standard periodic sets. Moreover, we allow degenerate (unbounded) polyhedra. The 
vertex-set of a polyhedron P, denoted by vert P, is defined as the set A n P. 

Theorem 3.1. I. Let V be a polyhedral subdivision ofW^ with a standard pe- 
riodic vertex-set A. Then the closure A('D) is a polyhedral cone in 5>q and 
A(D) is the set of all Q G 5>q satisfying 

(a) for every d-dimensional polyhedron P G T> the equalities 

{Nv,^, Q) = 0, 

for one (which can be chosen arbitrarily) affinely independent set ofd+\ 
vertices V C vert P and all w G vert P; 

(b) for every [d — l)-dimensional polyhedron F £ D the inequality 

Q) > 0, 

for one (which can be chosen arbitrarily) affinely independent set of d 
vertices V C vert F and two vertices w G vert P\F and w' G vert P' \ 
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F of the two adjacent d-dimensional polyhedra P,P' G T> with F = 

P n P'. 

II. The map T> i— > A(P) gives an isomorphism between the poset of Delone 
subdivisions of A. ordered by coarsening and the poset of closures of secondary 
cones ordered by inclusion. The closures of all secondary cones of Delone 
subdivisions form a polyhedral subdivision of S^q. 

Note that different choices of V, w and w' for the inequaUties could yield different 
conditions, respectively different forms iVyu{iu},ixi'- Nevertheless, these are the same 
on the linear subspace U defined by the equalities. In other words, their orthogonal 
projections Tru{Nvu{w},w') onto U are. all positive multiples of a uniquely determined 
form Nx)^F G U with {Nx>^f, N-d^f) = 1- 

Note also that for every G R'^ we have 

(4) Nv+v,w+v = Nv,w 

Therefore, and because the vertex-set of the subdivision is assumed to be periodic, the 
theorem gives only finitely many inequalities. This shows that A('D) is a polyhedral 
cone. 

Finally, let us remark that the theorem is valid for arbitrary periodic sets, that is, 
finite unions of lattice translates ti + L, if we replace 5>q by (^~^)*5>q(^~^) where 
A G GLrf(M) defines the lattice L = A'L'^. 

3.2. Proof of the fundamental theorem. In this section we prove Theorem 13. II We 
first give two propositions which both deal with the redundancies in the set of equations 
and inequalities, one would obtain by considering all possible choices of V, w and w'. 
Proposition 13.21 takes care of the equalities and Proposition [33]of the inequalities. The 
latter shows that the orthogonal projections of all the forms Ny^j^^^j ^^i for a facet F, 
onto the linear subspace defined by the equalities, are unique up to positive multiples. 

Proposition 3.2. Let Q £ S'^ and V C M"^ be an affinely independent set of cardinality 
d+1. Let c G M'^ and r > be such that Q[c — v] = for all v £ V. Then 

Q[w-c]-r^ = {Q,Nv,^). 

Proof. The proof is straightforward. We have 

Qlw-cl-r"^ = {Q,ww^) + {Q,-2wc*' + cc^) -r^ 

with a^, as in (|3]l. For each d G y we use the equality Q[v — c\ = which is 
equivalent to — (Q, vv^) = {Q, —2vc* + cc*) — r^. This yields the desired expression. 

□ 

Remark 3.3. By Proposition 13.21 the sign of {Q,Nv,w) has the following interpreta- 
tion: If it is positive, then w lies outside the circumsphere of the points in V, where the 
circumsphere is taken with respect to the norm induced by Q. If the sign is 0, then w 
lies on the circumsphere, and if it is negative, then w lies inside the circumsphere. In 
computational geometry this insph ere/outsph ere test is conveniently formulated using 
oriented mati-oid terminology (cf. iBVS+99l Chapter 1.8]): Let V = {vi, Vd+i) 
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be affinely independent points in with positive orientation and let w G W^. Then 
the chirotope 

■ 1 ... 1 1 



X(-ui,...,-Ud+i,-Lu)(Q) = sign 



Vi ... Vd+l w 

Q[vi] . . . Q[vd+i] Q[w] 



satisfies X{vi,...,Vd+i,n,)(.Q) = sign{Q , Nv,w) ■ 

Remark 3.4. Voronoi's theory and his description of secondary cones of Delone tri- 
angulations is based on linear forms Q[l,L') on S'^^ called regulators. Voronoi defines 
them for pairs of adjacent simplices {L,L') sharing a facet in a Delone triangulation. 
Let w' be the vertex of L' which is not a vertex of L. As in Proposition I3.2[ let V 
denote the vertex-set of L and define Ny^^'- Then Voronoi's regulator Q[l,l') is a 
positive multiple of {Ny^^i, ■). 

Proposition 3.5. Let P be a d-dimensional polyhedron in W^. Let U be the linear 
subspace of all Q ^ S'^ satisfying {Nv,w,Q) = for all affinely independent sets 
V C vert P of cardinality d+ \ and for all w G vert P. Let F be a facet of P. 

(1) Let V and V' be two sets of cardinality d containing affinely independent ver- 
tices of F and let w G vert P \ F and w' G W^. Then 

(2) Let V be a set of cardinality d containing affinely independent vertices of F 
and let u,w £ vert P \ F and w' G W^. Then 

Proof. In both cases we will show that the difference of the two considered forms lies 
in the orthogonal complement of U. 

(1) Every pair of affinely independent sets V,V' C vertF of cardinality d can 
be connected by a chain V = Vi, . . . ,Vn = V of affinely independent sets 
Vi C vert F of cardinality d such that | Vi n Fj+i | = d — 1. So we can assume 

n F'l = d — 1. Setting a.^,' = 1 there exist unique numbers and 
for u G y defining a affine dependency between the points w, w' and v in V. 
This defines the form A^yu{it)},u)'- 

We define vi by vi G V\V'. Since the affine hull of V and of V' equals the 
affine hull of F there exist numbers (3^ for v £ V' such that "^^v&v Pv = ^ 
and vi = J2v€V' f^f'"- Thus we have an affine dependency 

= aw + aw'+ ^ «ri + aviPv, 

v£V\{vi} veV 

= + a^'w' + y^ a^v + a^iPvV ' 

v£V\{vi} veV 

which defines the form iVy/u{ii,},io'- This gives 

Nvu{w},w' - Nv'U{w},w' = " 



Since ^ U y is a minimal affinely dependent set, we can choose an arbitrary 
vertex w of P and find that the right hand side is a multiple of A^yu{iu},iii- 
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(2) We take a closer look at the difference A'^^/ujiuj.ii;' — ^vu{u},w' and show it 
is a multiple of A^yu{iu},if Set a^' = 1 again and let and with v V 
be real numbers defining the affine dependency between the points w,w' and 
V in V. This defines the form Nyu^^y .^,,/. In the same way let a'^, = 1, a'^ 
and be real numbers defining Ny^j^u},w'- 

We set Pu = ck'^, /3w = —cnw and (3^ = a'^ — for v £ V. Then 
Pv'v = and J2 = where the sums run through all ?; G y U {u, w}. 
Thus 

(5) ^VU{w},w' — ^VU{u},w' = Pu^VU{w},u- 

n 

Note, if u and w' in the last calculation lie in opposite halfspaces with respect to 
the affine plane through V, then /?u > in Therefore, repeated application yields 
the following proposition, which we use for the proof of Theorem 13.11 and in Section 
l3.3l to prove Theorem 13.81 

Proposition 3.6. Let Vi , . . . , Vm C M.'^ be affinely independent sets of cardinality (i+ 1 
with \ Vi n = dfor z = 1, . . . , m — 1. Let w and Vi be on opposite sides of 
aff (Fj n ViJ^i)for i = 1, . . . , m — 1. Then 

m— 1 

Nv^,w = Nv^,'w + Yl (^i^V„-v,+i 
1=1 

with I'i+i G Vi+i \ Vi and positive constants Ui, for i = 1, . . . ,m — 1. 

With these propositions at hand, we can give a proof of the "fundamental theorem". 

Proof of Theorem \3J\ I. We show that A(D) is given by the set of listed linear equal- 
ities and inequalities. By (01) this implies that A(D) is a polyhedral cone, because the 
Delone subdivision induced by a Q G 5>q and a standai^d periodic vertex-set contains 
only finitely many Delone polyhedra up to Z*^ invariant translations. 

For Q G A('D) the linear equalities and inequalities are satisfied by Proposition 
EH 

Conversely, let us assume Q £ S'^ satisfies the linear equalities and inequalities for 
every polytope in V. By Proposition [3i2] and Proposition 13 . 5 1 we can assume that these 
are valid for all possible choices of V, v and w. Note that for the use of Proposition 
13. Sl it is crucial to observe, that for two linear- subspaces U, U' of S'^ with T = U OU' 
we have ttt = ttu o ttiji = tt^// o vrf/. 

Let P be a d-dimensional polyhedron in V. Let V C vert P be the vertex-set of a 
d-simplex. 

We show that all possible inequalities {Ny^^, Q) > 0, where w G P, are implied 
by the inequalities (A^v>i Q) ^ 0, where d G "D is either a vertex of P or a vertex of 
an adjacent Delone polyhedron of P. Assume w £ V\P. Then we choose a sequence 
of adjacent d-simplices with vertex-sets Vi, . . . , Vm in V satisfying the requirements 
of Proposition 13.61 In addition we require that each vertex-set is contained in the 
vertex-set of a fixed polyhedron of V, in particular Vi C vert P and w G Vm- Note 
that this can be achieved by looking at a refining triangulation of V, which we can 
choose arbitrarily. By our assumption on Q, we have {Nvi,v,+i , Q) > with equaUty 
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if and only if Vi, V^+i are subsets of the same polyiiedron. Thus by Proposition |3.6| and 
Proposition I3.2l we see that the equation {Ny^^w, Q) = imphes {Ny^^w, Q) > 0. 

Now we prove that a form Q ^ S'^ satisfying all linear conditions must be positive 
semidefinite. Let w and A G Z and write Xw = '}2veV ^^'^ ~ Z^dgV '^v- 
Then 

< {Nv,xn„Q) = A^QH - 

and hence Q [w] > since the coefficients depend only affinely on A. 

Furthermore, Q lies in the rational closure 5>q: The linear conditions are rational 
and thus define a polyhedral cone in which every extreme ray contains rational semi- 
definite quadratic forms. By the argument in the proof of Proposition 12.21 all quadratic 
forms in such a polyhedral cone lie in 5>q. 

We finally show that there exists a c G M'' and r > such that Q[w — c] > 
with equality if and only if w G vert P. A solution of the system of linear equations 
Q[v — c] = Q[v' — c], where v,v' V with V C vert P, gives c and r. Applying 
Proposition 13.21 shows that the inequalities Q[w — c] > ai^e valid with equality if 
and only if w G vert P. 

II. For the assertion on the isomorphism of posets given by the map V ^ •^(^)> 
we need to verify that the Delone subdivision I? is a true coarsening of V if and only 
if A('D) is strictly contained in A.(V'). This follows from I. because if P is a true 
coarsening of V we have equalities in the description of A(D) which are inequali- 
ties in the description of A.{V'). If on the other hand, A('D) is strictly contained in 
A (2?'), we know that all equalities and inequalities in the description of A(D') are 
also satisfied by elements of A(D) implying that P is a coarsening of V. Moreover, 
A(P) has to be contained in the boundary of A.{V'), because otherwise there would 
exist a Q G A('D') n A(D) implying V = V'. Thus at least one of the inequalities in 
the description of A(D') is fulfilled with equality for the elements of A(D). Hence V 
is a true coarsening of V. 

Since every form in 5>q defines a unique Delone subdivision of A, 5>q is subdi- 
vided into secondary cones. This subdivision has the face-to-face property, hence is a 
polyhedral subdivision, due to the poset isomorphism: Otherwise we would find De- 
lone subdivisions V and V, such that A(D) and A.(V') have the same dimension and 
such that there is a common relative interior point, but with A(P) / A(D'). Thus Q 
would yield two different Delone subdivisions, which is not possible. □ 

3.3. Flips and bistellar neighbors. Given a linear subspace T of S'^, we know by 
Theorem 13.11 that 5>q H T is covered by the topological closures of T-generic T- 
secondary cones. As before we assume that the vertex-set of the considered Delone 
subdivisions is a fixed standard periodic set A. Given a T-generic Delone subdivision 
V and its T-secondary cone A.TiT>) we want to determine its contiguous T-secondary 
cones. For this let .7^ be a facet of Ar('D) which is not a dead-end. Let V' be the 
T-generic Delone subdivision with A(P') being contiguous to A('D) at J^, hence 
with J" = At(P) n At(P')- The transition from V to V is called a T-flip; the 
subdivisions V and V are refeiTcd to as bistellar neighbors. 

In order to describe the T-flip, let Njr G S'^ denote the form which is uniquely 
determined by the following conditions: 
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(i) NjreT with {Njr, N^f) = 1, 

(ii) (iV^, g) = for all Q G 7", 

(ill) {N^, Q) > for all Q G At{V), 

(iv) {N^, Q) <0 for all Q G At(P')- 
We collect facets F of P whose forms Nx>^f (as defined after Theorem |3.1| ) are positive 
multiples of Njr when projected onto T: 

njr = {F eV : dimF = d - 1 and TVjr = a • ttt{Nv,f) for a a > 0}. 
On TZjr v/e define an equivalence relation by 

{there exist Fi, . . . , F„ G 7^J^ with F = Fi, F' = F^, 
and d-dimensional polyhedra Pq, . . . , Pn G V 
with Fi = Pj_i n Pi for i = 1, ... ,n. 

Thus by definition, the facets in each equivalence class are connected by a chain of 
adjacent d-dimensional polyhedra in V. The union of these polyhedra is a polyhedron 
again: 

Proposition 3.7. Let D be a Delone subdivision and C an equivalence class of (O. 
Let Vc be the set of all vertices of d-dimensional polytopes in T) with a facet in C. Then 
(convVc) n A = Vc- 

A proof of the proposition is given at the end of this section. It, as well as the 
description of T-flips can conveniently be given using the lifting map wq : M*^ — > 
M'^+^ defined by wq{x) = {x,Q[x]) for a Q G 5>q. The polyhedra convVc in 
Proposition 13.71 are called repartitioning polyhedra of T. For a Q G Ar('D) and a 
repartitioning polyhedron P we define the [d + 1) -dimensional polyhedron 

P{Q) = conv{wQ{v) : v G vertP}. 

For each of its facets F we have an outer normal vector n{F) G M"'"'''^. The facets 
are divided into three groups. We speak of a lower facet, if the last coordinate of the 
normal vector satisfies {n{F))ii+i < 0; we speak of an upper facet, if {n{F))d+i > 
and of a lateral facet, if {n{F))d+i = 0. We show below that these notions are 
independent of the particular choice of Q G A.t{'D). 

For the description of the T-flip in the following theorem, let tt : M'^"'"^ ^ 
denote the projection onto the first d coordinates. Figure 1. gives an example for the 
change from "lower to upper hull" in the case of a planar repartitioning polytope. Note 
that the situation in higher dimensions can be much more complicated than the picture 
might suggest. 
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Figure 1 . Lifting interpretation of "upper to lower hull change". 

Theorem 3.8. Let T be a linear subspace of S'^ and J- be a common facet of two 
contiguous T -secondary cones of T -generic Delone subdivisions T> and T>'. Then we 
obtain the d-dimensional polytopes of D' from those ofD by choosing an arbitrary 
Q G and by doing the following for each repartitioning polytope P of T: 

(i) We remove all polytopes it{F) of lower facets F of P{Q). 

(ii) We add all polytopes tt{F) of upper facets F of P{Q). 

Remark 3.9. Our theory is similar to the theory of equivariant secondary polytopes of 
regular polytopal subdivisions, as recently described by Reiner in IIRei02ll . With slight 
modifications our explicit description of the flipping procedure works in his setting as 
well. The linear subspace T is in this case containing all PQFs which are invariant 
under a given finite subgroup of GL(i(Z) (see Section |4|. 

As a preparation for the proofs of Theorem 1 3 . 8 1 and Proposition l3.7[ consider the set 
(7) conv{wQ{x) : X e A}. 

It is a locally finite polyhedron, meaning that the intersection with a polytope (bounded 
polyhedron) is a polytope again. A first and important observation is that the facets of 
the polyhedron ^ yield the li-polytopes of the Delone subdivision when projected by 
vr onto the first d coordinates (see Figure 2): 

Proposition 3.10. Let Abe a standard periodic set and Q G S^q- Then 

Del{Q) = {vr(F) : F face of ^ of dimension less than or equal d}. 




Figure 2. Delone subdivision obtained from lifting. 

Proof. For a polytope P ^ V and a set y C vert P of cardinality d + 1 and with 
affinely independent vertices, the set of all {x,y) G M'^+^ with 

(8) y = Q[x]- {Nv,^, Q) = Y1 «-^[^] 

vev 

is an affine hyperplane in R'^+^ . Recall from the definition of Ny^x in © that the a-^ are 
uniquely defined by the affine dependency x = J2vgV ^"""^ ^ith ^^gy/ ~ ^- Thus 
we see that the hyperplane given by ([8]l contains the d+1 lifted points wq{v). More- 
over, it is a supporting hyperplane of the polyhedron ([7]) if and only if {Nv,x, Q) > 



14 



MATHIEU DUTOUR SIKIRIC, ACHILL SCHURMANN, AND FRANK VALLENTIN 



for all a; G A. By Proposition 13.21 this is the case if and only if Q[x — c] > for a 
suitable c G M'^ and r > 0. Here, equality holds if and only if a; G vert P. □ 

The hyperplane given by equation ([8]l in the proof is a different one for subsets V 
and V of different Delone polytopes. The forms Nxi^f relate the change for adjacent 
d-poly topes P and P' in V with common facet F = PnP'. If the segment connecting 
X with a vertex v of P intersects facets Fi, Fn of V, then by Proposition 13.61 and 
by the definition of N-p^p (after and because of Theorem 13. II ) we have 

n 

(9) y = Q[x] - {Nv,^, Q) = Q[x] - Pi{Nv,Fo Q) 

1=1 

with suitable positive constants depending on x. 

Equation ^ yields simple proofs for Proposition 13 .7 1 and Theorem [331 

Proof of Proposition \3. 71 Due to ^ and the definition of C, all points wq{v) with 
G Vc are coplanar in M'^^^. Suppose the union of the d-dimensional polyhedra of 
V with a facet in C is strictly contained in conv Vc- Then there exists a facet F of one 
of the polyhedra such that F ^ C and F does not belong to the boundary of conv Vc- 
Hence, there exist two vertices w,w' G Vc in opposite half spaces with respect to 
aff F. Then by Q and since {Nx)^f, Q) > we see that wq{w), wq{w') and wq{v) 
with 17 G F can not be coplanar, which is a contradiction. □ 

Proof of Theorem]!^ Choose Q G A(P) and Q' G ^{V) and let P be a repartition- 
ing polytope of = A(P)n A(D'). Then by formula © and since Nv,f = -Nv',f 
for all facets F contained in P, the upper facets of P{Q) project to the same polyhe- 
dral subdivision of P as the lower facets of P{Q') and vice versa. The lower facets 
give the corresponding Delone subdivisions by Proposition 13.101 which proves the as- 
sertion. □ 

4. FiNITENESS AND EQUIVARIANCE 

For the sake of simplicity and because it is sufficient for our applications in Section 
[Hand Section [9] we restrict ourselves from now on to the case of Delone subdivisions 
with vertex-set U^. A more general discussion can be found in IISch07l . In particular it 
is shown that the results of this Section extend to the case of rational standard periodic 
vertex-sets A = U™ ^ U + with G Q*^ for i = 1, . . . , m. 

Let T be a linear- subspace of S'^. We say that two T-secondary cones At and 
are T-equivalent if there is a <^ G GLd(Z) so that g^A^g = A^ and g^Tg = T. 
Otherwise we say that they are T-inequivalent. In this section we discuss assumptions 
on T which ensure that there exist only finitely many, T-inequivalent T-secondary 
cones. Note that Voronoi's classical theory (cf. Theorem 12.31 ) deals with the case 
T = 5'^. 

The main difference with the classical theory is that for general T there may be 
infinitely many, T-inequivalent T-secondaiy cones: 

Example 4.1. Consider a rational PQF Q with trivial automorphism group. The orbit 
of Q under the action of GLd{1j) contains only rational PQFs. We choose a subspace 
T through Q in which the only rational forms are multiples of Q. So T does not 
contain another element from the orbit of Q. Thus the setwise stabilizer of T is trivial. 
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In this situation there are infinitely many T-generic secondary cones whenever the set 
Tn5>Q is not closed in S"^. Then the intersection of T with the boundary of 5>q is not 
covered by finitely many T-dead-ends. The maximum dimension of the intersection 
of S^Q with a supporting hyperplane of S^^ is C^^^) - d (see IIBR79I §10]). Thus 
the intersection of T with the boundary of 5>q cannot be covered by T-dead-ends if 

dim(T n 5>o) > ("'t^) —d + 1. In particular for all li > 3 there exist examples of 
subspaces T with infinitely many, inequivalent T-generic cones. 

In the following we show that there are only finitely many, T-inequivalent T-secondary 
cones if T is the linear subspace which is stabilized pointwise by a finite subgroup G 
of GLd(Z). We develop this equivariant theory in analogy to the theory of G-perfect 
forms of Berge, Martinet and Sigrist IIBMS92II . In Remark l4!4l we give a unifying view 
on both theories. 

First let us recall some definitions and basic results (see IIBNZ73I Section 2]). Let 
G be a finite subgroup of GL(i(Z). The linear subspace 

J^{G) = {Q G 5^ : g^Qg = Q for all g e G} 

is called the space of invariant forms of G. The pointwise stabilizer of J^{G), which 
we denote by 

B{G) = {g& GLrf(Z) : g'Qg = Q for all Q G T{G)}, 

is called the Bravais group of G. Note that B{G) can be strictly larger than G. The nor- 
malizer of a subgroup G of is defined by N{G) = {n S GL£^(Z) : n~^Gn = 

G}. One important property of N{B{G)) is that it is the setwise stabilizer of the sub- 
space of invariant forms 

N{B{G)) = {<? G GL,(Z) : g'HG)g = T{G)}. 

A proof can be found for example in ]Jaq95[ Lemme 3.2]. 

Now we are ready to state the main result of this section. For this recall that an 
open polyhedral subdivision is a non-intersecting decomposition into polyhedral cones 
which are open with respect to their affine hull. 

Theorem 4.2, Let V he an open polyhedral suhdivision of S^q on which GL^i'^) acts 
by {g, A) g^Ag. Suppose that this action gives only finitely many orbits. Define a 
polyhedral subdivision of !F{G) n by 

VriG) = {A n HG) : A G P}. 

TJien the normalizer N{B{G)) acts on 'Pjf(g') and there exist only finitely many orbits 
with respect to this action. 

This together with Theorem |23] gives the following corollary, which completes our 
equivariant version of Voronoi's theory. 

Corollary 4.3. Let G be a finite subgroup o/GLrf(Z) cmd let T = J-{G) be the space 
of invariant forms. Then there exist only finitely many T-inequivalent T-secondary 
cones of D clone subdivisions ofU^. 

Proof of Theorem \4~2\ We start with a definition. For A G P we define the automor- 
phism group Aut(A) = {g e GU(Z) : g^Ag = A}. 
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We consider the set 

Oa,g = {g'^g n HG) ■■ g G GLd(Z) and g'^g n HG) + 0}. 

The normahzer N{B{G)) stabihzes J^{G) setwise, and hence it acts on Oa,g- We 
show that Oa,g is a finite union of A^(;B(G)) -orbits, that is, there are gi, . . . , gm G 
GLd(Z)with ' 

m 

Oa,g = [jWigjAgi n TiG))n : n G N{B{G))}. 

i=l 

Then the statement of the theorem follows because the set {Oa,g '• A G V} is finite 
by the assumption made on V. 

Let gi,g2 be in GL(i(Z) so that gjAgi n J^{G) G Oa,g for z = 1, 2. The group 
B{G) is a subgroup of Aut{g\Agi) and Aut((72A52) which can be seen as follows. 
By definition the group B{G) stabilizes the set g\A.gir\lF{G) pointwise. Since GLrf(Z) 
operates on V and gjAg-i n T{G) ^ every element of B{G) has to stabilize gjAg-i 
setwise. 

Hence, giB{G)gi^ and g2B{G)g2^ ai^e subgroups of Aut(A). Assume that 
g2B{G)g2^ ■ Then n = g2^gi is an element of the normalizer N{B{G)). Furthermore, 
we have 

n*(ff* A<72 n J^{G))n = n*<7* Agsn n n'j^{G)n = g{Agi n J^{G) 

because N{B{G)) stabilizes T{G) setwise. Thus different N{B{G)) -orbits in Oa,g 
induce different subgroups of Aut(A). 

For every A £ V the group Aut(A) is finite (This fact was already proved by 
Nakamura in IINak75[ Lemma 1.2]. For completeness we give an argument, which also 
is a bit more elementary than Nakamura's.): Let TZi, . . . ,71^ G 5>q be the extreme 
rays spanning the closed polyhedral cone A. Choose forms Qi G TZi in the following 
way: If for a pair TZj there is a [/ G Aut(A) with U^TZiU = TZj then U^QiU = 
Qj. We have A = cone{Qi, . . . , Qk} and Q = Yli=i Qi G A and in particular 
Q G 5>o- So, Aut((5) = {g £ GU(Z) : g^Qg = Q} is finite. By construction 
Aut((5) contains Aut(A). So, the group Aut(A) has only a finite number of different 
subgroups. So, Oa,g is a finite union of A^(i3(G))-orbits. □ 

We close this section with two remarks. 



Remark 4.4. The proof of Theorem 14.21 is essentially an adaptation of proofs of 
Jaquet-Chiffelle ( |Jaq95[ Theoreme 5.2]) and Opgenorth (| |Opg95[ Theorem 4.2.18]). 



They show that there exist only finitely many inequivalent G-perfect forms (definition 
is given below). Proofs for their theorems can be derived from Theorem 14.21 in the 
following way: 

Let m be a positive number and let Vm be the set 

Vm. = {Q G 5>o : > m for all v€Z'^\ {0}}. 

This is a convex, locally finite polyhedral cone. For each face F of Vm we define 
the (relatively) open polyhedral cone Ap = {XQ : A > and Q G relintF}. Then 
V = {Ap : F face of Vm} gives an open polyhedral subdivision of S^q as required 
in Theorem 14.21 
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If G is a finite subgroup of GLrf(Z), tiien Q is called G -perfect if for the cone A.p 
with Q G Ai? we have dim(Ai? n J^{G)) = 1. Two G-perfect forms Q and Q' 
are called G -equivalent if there is a (7 € N{B{G)) and a positive A so that g^Qg = 
XQ'. Otherwise they are called G-inequivalent. By Voronoi's first memoir IIVor07l the 
polyhedral subdivision V satisfies the assertion of Theorem 14.21 Hence, there are only 
finitely many G-inequivalent G-perfect forms by our Theorem l4.2l 

Usually G-perfect forms (cf. IIMar03[ IBMS921 |Jaq95| |Opg95| |Qpg01[ ) are defined 



via normal cones of faces Fq of Vm n !F{G) in T{G). A face F of Vm is uniquely 
characterized by the set 

Min(F) = {v ^ if : Q{v] = m for all Q S F}, 

which is independent of the particular choice of m. The normal cone M{F) of F 
is given by Af{F) = cone{vv'^ : v G Min(F)} and the normal cone of the face 
Fg = F n T{G) in T{G) is obtained by an orthogonal projection of J\f{F) onto 
T{G). If different inner products are used, the resulting cones differ as seen in the 
cases of |Jaq95[ and |Qpg95| |. 



Remark 4.5. The case when the linear subspace T is pointwise stabilized by a fi- 
nite subgroup of GL(i(Z) is not the only case when there are only finitely many T- 
inequivalent T-secondary cones. Bayer-Fluckiger and Nebe ( IIBN05I Theorem 3.1]) 
give another sufficient condition for T which ensures finiteness in the two-dimensional 
case. However, we are not aware of further results regarding this question. For all the 
subspaces T we considered for the applications, there were only finitely many T- 
inequivalent T-secondary cones. These subspaces are all spanned by rational forms, 
but this alone does not guarantee finiteness for > 3, as shown by the following 
example due to Yves Benoist (private communication): 

Consider the subspace T spanned by + 2y^ + and xy. It contains the two 
positive semidefinite, but non-positive forms (x ± \piyf' + . In any small neigh- 
borhood of these two forms, we find infinitely many forms in T n S'^^ with pairwise 
differing Delone subdivisions. In particular, the minimal non-zero vectors v ^ 
tend towards the kernel of one of the two forms; the Delone subdivisions have dif- 
fering edges convjO,!)} with v = {p,q,OY and | sufficiently close to =F\/2- Thus 
there are infinitely many different T-secondary cones, but the setwise stabilizer of T 
in GLrf(Z) is finite. In order to see this, note that the three lines of rank-2 forms in T 
have to be permuted. Note that the same argument does not apply to dimension 2. The 
space spanned by x^ + 2y^ and xy has an infinite stabilizer. 

5. Other generalizations of Voronoi's Reduction Theory 

In the seventies, Mumford and Satake described a general procedure for compacti- 
fying quotients of bounded symmetric domains by arithmetic groups (cf. MAMRTTSID . 
Their constructions use special decompositions of the domains of interest into ratio- 
nal polyhedra. As an application of this general method, Namikawa ( IINam76l ) con- 
structed a compactification of Siegel modular varieties using Voronoi's reduction the- 
ory. Siegel modular varieties ai^e moduli spaces of Abelian varieties. Over the complex 
numbers they arise as quotient of the Siegel upper halfspace by an arithmetic group. 
Namikawa's construction opened the possibility to understand the boundary of the 
moduli space with help of degenerations. Later algebraic geometers refined this work 
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and they used also generalizations of Voronoi's reduction theory for this. In this section 
we want to review these generalizations and compare them to ours. For more informa- 
tion about the geometry of Siegel modular varieties we refer to the survey ||HS02i of 
Hulek and Sankaran. 

Oda and Seshadri study in IIQS79I a generalization of Delone subdivisions which 
they call Namikawa decompositions: Let be a Euclidean space where an orthogonal 
decomposition into subspaces E' and E" is given: E = E' 1. E", and let A C 
be a lattice so that A n £" is also a lattice. For a vector t/j G E" the Namikawa 
decomposition is 

D^ = {7r'{D{x)):xeE' + ^}, 
where vr' :£'—>£" is the orthogonal projection onto E' and 

D{x) = conv{^i, . . .,^r}, 

with vertices G A closest to x, is a Delone polytope of A. Then, Oda and Seshadri 
explain in IIOS791 Proposition 2.3] how the Namikawa decomposition changes when 
one varies the vector ip along E". Our generalization of Voronoi's reduction theory 
goes into a different direction: We fix a lattice A in a real vector space of finite dimen- 
sion and a basis of A and study how the Delone subdivision of A changes when one 
varies the inner product along a given subspace (in the space of Gram matrices). 

Alexeev and collaborators (cf. IIABH021 lAle02l IAle04ll and the survey IIAle06[ Sec- 
tion 5]) consider so-called semi-Delaunay decompositions: Let M C M'^ be a lattice 
and r C M be a sublattice of finite index. A F-periodic polyhedral subdivision of 
is called a semi-Delaunay decomposition when it is the projection of the lower hull 
(see Section [331 ) of the lifted points (m, /i(m)) with m £ M and /i : M — > M a 
function of the form h{m) = q{m) + r{m) where g : M'^ ^ M is a positive semidef- 
inite quadratic form and r : M/T — > M an arbitrary function on the cosets. Delone 
subdivisions of rational standard periodic sets, i.e. sets of the form ([Hi with rational 
t[, are covered by this construction: Take T = 1/ and an overlattice M containing all 
t[. Define r(i'- + F) = and r(m + F) large enough for all m & M not in the rational 
standard periodic set. 

Nevertheless, we believe that our treatment has its merits since it is very explicit and 
can be immediately used for computations. Moreover, our generalization of Voronoi's 
reduction theory to positive definite quadratic forms with prescribed automorphism 
group, which is essential for our applications, is not covered by the reviewed construc- 
tions. 

6. Optimizing lattice sphere coverings 

The lattice covering problem is a classical problem in "Geometry of Numbers". 
Roughly speaking, the problem is concerned with the determination of the most eco- 
nomical way to cover W^. We deal with it in Section [8] and [9] where we apply our 
generalization of Voronoi's theory. In this section we give the necessary definitions 
and describe briefly the methods used for the applications. For further reading and 
more background we refer the interested reader to IISV06II . 

6.1. Definitions. We assume that L is a lattice of full rank d, that is, there exists a 
matrix B G GLd(IR) with L = BZ"^. The determinant det(L) = | det{B)\ > of L is 
well defined and does not depend on the chosen basis. Let B'^ = {a; G M'^ : ||a;|| < 1} 
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denote the solid unit sphere in M°'. The Minkowski sum L + aB'^ = {v + ax : v G 
L,x ^ B"^}, with a G ]R>o> is called a lattice covering if M"^ = L + aB'^. The 
covering radius fi{L) of L is given by 

= min{/i : L + fiB'^ is a lattice covering}. 

For a £ M we have 

fi{aL) = |a|/i(L) 

Thus, the covering density 

with Kd = vol i?*^, is invariant with respect to scaling of L. Note also that is an 
invariant of the isometry classes. 

The lattice covering problem asks to minimize among all d-dimensional lattices. 
It has been solved only for dimensions d < 5 based on Voronoi's classical reduction 
theory and the knowledge of all Delone subdivisions in these dimensions. 

If we work with PQFs, the definitions above translate in the following way: 



Here, the inhomogeneous minimum fi{Q) is given by 

fj.{Q) = max mill Q[x — v], 

satisfying det(L) = y^det{Q), fi{L) = y/fiiQ), for a corresponding lattice L ob- 
tained from Q. 

Given a linear subspace T C S'^, let 

Gt = inf e(Q) 

denote the bound on the covering density with respect to T. This bounds is computed 
in our applications described in Section[8]and Section|9] 

6.2. Determinant maximization problems. Let T C cS'' be a linear subspace. To 
compute @T, we consider all inequivalent T-generic Delone subdivision V (if possi- 
ble) and minimize G among all PQFs in Ar(D). This can be achieved by solving a 
determinant maximization problem. These ai^e convex programming problems of the 
form 

minimize c*a; — logdetG(a;) 
^ ^ subject to G{x) >~ 0, F{x) >z 0, 

over the variable vector x G M^. The objective function contains a linear part given 
by c G M^. The affine maps G : ^ M™><™, as well as F : ^ M"^", are 
given by 

G{x) = Go + xiGi + --- + xdGd, 
F{x) = Fo + xiFi + --- + xdFd, 

where Gi G M™^™ and Fi G M"^", i = 0, ... ,D, are symmetric matrices. The 
notation G{x) >- and F{x) ^ gives the constraints ''G{x) is positive definite" 
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and "F{x) is positive semidefinite". Note that we specialize to a semidefinite pro- 
gramming problem, if G{x) is the identity matrix for all x G MP . One nice feature of 
determinant maximization problems is that there is a duality theory similar to the one 
of linear programming (see IIVBW98II ). It allows us to compute an interval in which 
the optimum is attained, the so called duality gap. 

We can express the condition /u(Q) < 1 as a linear matrix inequality (LMI) F{x) ^ 
where the optimization vector x is given by the coefficients of Q = Yl,i=i ^i^i' 
with respect to some basis {Ai , . . . , ^dim t) of T. This is seen by the following propo- 
sition due to Delone et al IIDDRS70I (cf. IISV06I Proposition 7.1]), together with the 
crucial observation that the inner product (•, •) defined by {y, z) = y^Qz can be ex- 
pressed as 



{y,z) = {Q,yz') 



dimT 

E 

i=l 



Xi{Ai,yz^ 





/ 4 




Vi) 




■ {Vd,Vd)\ 




{vi,vi) 




Vi) 


{Vl,V2) . 


■ {Vl,Vd) 


BRl(Q) = 


{V2,V2) 




Vi) 


{V2,V2) . 


■ {V2,Vd) 






{Vd 


Vi) 


{Vd,V2) . 





hence as a linear combination of the parameters Xj. 

Proposition 6.1. Let L = conv{0, vi, . . . , Vd^ Q be a d-dimensional simplex. 
Then the circumradius of L is at most 1 with respect to (•, •) if and only if 



>- 0. 



Since a block matrix is semidefinite if and only if the blocks are semidefinite, we 
get the desired result. Here, we call two Delone polytopes P and P' equivalent with 
respect to Q ^ '5>0' if there exists a [/ in the automorphism group 

Aut(Q) = {U e GLd(Z) : U^QU = Q] 

of Q and a G Z'^ with P = v + UP'. 

Proposition 6.2. Let Q G 5>g be a PQF with Delone subdivision V. Let Pi , . . . , P„ 

be a representative system of d-dimensional Delone polytopes in T>, which are in- 
equivalent with respect to Q. For every Pi choose a d-dimensional simplex Li with 
vert Li C vert Pi. Then 

\ 



li{Q) < 1 



BRlAQ) 











BRl,(Q) 













v 















^ 0. 



Thus the constraint fi{Q) < 1 can be brought into one LMI of type F{x) ^ 0. 
We can add linear constraints on the parameters Xj by extending F by a 1 x 1 block 
matrix for each linear inequality. In this way we obtain one LMI for the two constraints 
fJ'iQ) < 1 and Q G A(P). Therefore we can determine a PQF in Ar(X>) minimizing 
the covering density by solving the following determinant maximization problem. 

minimize — log det ( Q ) 

subject to Q G Ar(P), niQ) < 1. 
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Here F{x) is as described above, and c = 0, and G{x) = X^fl™^ XjAj in (ITOl ) where 
the Ai form a basis of T. 

Remark 6.3. The problem is a convex programming problem. Hence (see e.g. IISV061 
Proposition 9.1]), if the automorphism group of the considered Delone subdivision V 
contains the pointwise stabilizer of the subspace T, we know that a PQF minimizing 
the covering density in A.(V) is contained in T . Hence by choosing the subspace T 
carefully we sometimes can dramatically decrease the dimension of the optimization 
problem. 

Remark 6.4. Determinant maximization problems have not only theoretically nice 
features, but also can be solved in practice quite well, for example with interior point 
methods. For the solution of the optimization problems described above we use the 
software package MAXDET irWVB96ll of Wu, Vandenberghe, and Boyd as a subroutine. 
By using rational approximations it is possible to give mathematical rigorous error 
bounds or even proofs of local optimality. For details we refer to IISV06II . Utilizing 
these ideas we developed a program rmd (rigorous maxdet) and based on it coop, 
a covering optimizer. These programs, together with a short tutorial can be obtained 
from our web page BSVOSall . They allow to approximate — mathematically rigorously 
— the optimum covering density with respect to a given T-secondaiy cone. 

7. Algorithmic issues 

In this section we discuss several algorithmic issues which have to be resolved to 
turn our generalization of Voronoi's reduction theory into an effective procedure. As 
in Section |4] we assume that all Delone subdivisions have as a vertex-set. For the 
original theory of Voronoi a similar discussion can be found in IISV061 Section 5.3]. 
Here we emphasize those points where the generalization differs from the classical 
case. First we adapt IISV06I Algorithm 1], which enumerates all inequivalent Delone 
triangulations, to an algorithm enumerating all T-inequivalent T-generic Delone sub- 
divisions. As in Example 14. 1 [ our Algorithm 1 below does not necessarily stop after 
finitely many steps, depending on T. 

Apart from this, when we compaie Algorithm 1 to the one of the original theory, 
essentially two new tasks arise: 

(1) We have to find an initial T-generic Delone subdivision (cf. Algorithm 2), 
whereas in the classical case this can be given by "Voronoi's first subdivision". 

(2) We have to check whether two T-generic Delone subdivisions are T-equivalent 
(cf. Algorithm 3). 

For Algorithm 2 we need one more definition: The dimension of the linear span 
of At('D) is called T -rigidity index. For the classical setting T = S'^ the rigidity 
index was introduced by Baranovskii and Grishukhin in MBGOH . Algorithm 2 is a 
randomized algorithm and belongs to the class of so-called Las Vegas algorithms. That 
is, it always produces coiTcct results whereas the running time is a random vaiiable. 

Algorithm 3 checks whether two T-generic T-secondary cones are T-equivalent. 
For this we need a definition. 

Definition 7.1. We say that a positive semidefinite quadratic form Q is rational nor- 
malized if it is integral and the entries have greatest common divisor equal to one. 
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Input: Lineal- subspace T Q S'^ 

Output: Set IZ of all T-inequivalent T-generic Delone subdivisions 
Y <— {Vi}, where Vi is a T-generic Delone subdivision (Algorithm 2). 

while there is a 2? G y do 

Y ^ Y\{v}. n^nu {V}. 

Compute the Unear inequalities of Ay (2?) as in Theorem [3TT] 
Compute the facets Fi, . . . , F„ of At('D). 
for i = 1, . . . , n do 
if Fi is not a dead-end then 
Compute the bistellar neighbor Vi of V, defined by Fi as in Theorem l3.8l 
if Vi is not T-equivalent to a Delone subdivision in the set 
n U {Vj : j e Fj not a dead-end} (Algorithm 3) then 

Y^YU{V,}. 
end if 
end if 
end for 
end while 

Algorithm 1. Enumeration of all T-inequivalent T-generic Delone subdivisions. 



Input: Linear subspace T C S'^ 

Output: A T-generic Delone subdivision Vi 

repeat 

Choose a random Q in T n 5^q. 
Compute Vi = Del(Q) (see IIDSV07II ). 

Compute the T-rigidity index m of Vi by the linear equalities in Theorem |3TTJ 
until m = dim T. 

Algorithm 2. Finding a T-generic Delone subdivision. 

For Algorithm 3 we have to test whether two PQFs are equivalent and we have to 
compute the automorphism group of a given PQF. For both tasks exist efficient algo- 
rithms by Plesken and Souvignier MPS 971 . Implementations are part of the computa- 
tional algebra system MAGMA and of the computer package CARAT (Crystallographic 
AlgoRithms And Tables, cf. IIOPS98II ). 

The closed polyhedral cones At{V) are rational (assuming that the vertex-set is 
Z*^). That is, they have a description by rational inequalities (cf. Section [3]) as well 
as a description A.t{V) = cone{i?i, . . . , Rk} as a cone generated by finitely many 
rational normalized quadratic forms lying in extreme rays of A.t{V). 

Dealing with symmetries of a polyhedral cone, the characteristic form Q = X^JLj^ Ri 
is of great importance. This is due to the following proposition which is used in Al- 
gorithm 3. We call two rational polyhedral cones V and V' in S'^ equivalent, if there 
exists aU G GldC^) such that U^VU = V. As in the proof of Theorem |42l the 
automorphism group of V is defined by Aut('P) = {[/ G GL(i(Z) : U^VU = V}. 
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Input: Lineal" subspace T C S'^, and T-generic Delone subdivisions Vi, V2 
Output: Yes, if Pi and D2 are T-equivalent; No, otherwise 

TZi ^ {rational normalized element R of extreme rays of AxiVi)}, for i = 1,2. 
Qi EfleTe, RJori = 1,2. 
Compute Aut((5i). 

if there exists A £ GU(Z) and B £ Aut((5i) 

with A^QiA = Q2 and {BAyT{BA) = T then 

return Yes. 
else 

return No. 

end if 

Algorithm 3. Checking T-equivalence. 



Lemma 7.2. The automorphism group of a characteristic form contains the automor- 
phism group of its polyhedral cone. The characteristic forms of two equivalent rational 
polyhedral cones are equivalent. 

Proof. Suppose that for two polyhedral cones P, -P' C there exists a [/ G GLrf(Z) 
such that U*VU = V' . Extreme rays of V are mapped onto extreme rays of V' . We 
need to show that this is also true for the uniquely determined rational normalized 
elements of the extreme rays. For this let Q be such a rational normalized element of 
an extreme ray of V. It is mapped onto the integral matrix U^QU = aQ' defining an 
extreme ray of V' , with rational normalized Q' . Thus a G N. On the other hand Q' is 
mapped onto an integral multiple of Q via U^^ and (C/~'^)*(5'(f/~^) = (l/a)Q. Thus 
a = 1, which proves the assertion. □ 



8. Application L Classification of totally real thin number fields 
In this section we apply our theory to a problem in algebraic number theory. Re- 



cently, Bayer-Fluckiger introduced in |Bay06| the notion of thin algebraic number 



fields. A thin algebraic number field is Euclidean due to a special geometric reason, 
which we explain in Section 18.11 She proved that there exist only finitely many thin 
number fields (| Bay06[ Proposition 11.4]). 



In IIBN05II Bayer-Fluckiger and Nebe gave a complete list ( IIBN051 Theorem 5.1]) 
of 17 candidates for totally real thin algebraic number fields. They proved that 13 of 
them are thin and one of them is weakly thin but not thin. However, in three cases 
they did not know how to decide whether the fields are thin or not. Using our theory 
we finish the classification of totally real thin algebraic number fields. In particular we 
show in Section [8^ that the three open cases do not give weakly thin fields. 



8.1. Background and definitions. Let us recall some standard definitions from alge- 
braic number theory. 

Let K be a number field of degree n = [K : Q]. From now on we assume that K 
is totally real, that is, for all embeddings of fields ai : K —>■ C with i = 1, . . . , n we 
have ai{K) C M. 
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Let Ok be its ring of integers which is a Z-module of rank n. Thus we can write 
Ok = Z'j-'i + • • • + '^LOn- We embed Ok into M" via the map a defined by a{x) = 
{ai{x), . . . , cr„(x)). Hence, (j{ok) is a lattice of rank n. 

By V we denote the set of a e K with <Tj(a) > for alH = 1, . . . , n. Then 

n n 

{x,y)a = ^(7i{axy) = '^ai{a)ai{x)ai{y) 
1=1 1=1 

defines an inner product on K and so on M". Every inner product (•, •) : K x — > M 
with the additional property {x, yz) = {xz, y) for all x, y, z G K is of this form. 

We denote by KK,a the lattice which is given by the pair {ok, {■, •)«) where a ^ V. 
We say that K is weakly thin if for @{K) = mina^-p @{AK,a) we have the inequality 



/ n" 

and we say that K is thin if we have strict inequality. By t{K) we denote the number 
on the right hand side of this inequality. 

Let us briefly explain the motivation of this definition. First note that (det A^^ i)^ 
equals the discriminant cIk of K. Let : — > M be the norm of K which is given by 
= Yli=i ^i{x)- A number field K is called Euclidean if its ring of integers Ok is 
an Euclidean ring with respect to the absolute value of the norm function. Equivalently, 
K is Euclidean if and only if its Euclidean minimum 

M{K) = sup inf \N{x-y)\ 
x(zKy<^°K 

is strictly less than 1. Currently there is no algorithm known which computes the 
Euclidean minimum of a number field. It is also an open problem if there exists a 
finite or an infinite number of Euclidean number fields. In IILen77ll Lenstra, motivated 
by work of Hurwitz, gave several bounds on the Euclidean minimum of a number field 



using methods from "Geometry of Numbers". In |Bay06| Bayer-Fluckiger extended 
these results. By using the inequality between the arithmetic and geometric means to 
relate the norm function N to the norm ||x||a = {x, x)a = ■\/Yll=i ^iC^^^^), she 



showed that thin fields are Euclidean fields (| Bay06[ Proposition 1 1.2]) and that there 



exist only finitely many thin fields (| Bay06[ Proposition 1 1.4]) 



8.2. Classification. By applying the lower bound for sphere coverings of Coxeter, 
Few and Rogers IICFR59II . Bayer-Fluckiger showed that the degree of a totally real 
thin number field is at most 5 (| Bay06[ Proof of Proposition 11.4]). All number fields 



having low degree and low discriminant are known ( IIPH05I ). Using this list Bayer- 
Fluckiger and Nebe gave a complete list of 17 candidates of totally real number fields 
which might be thin. In Table 1 we list these candidates K together with the relevant 
parameters: The degree n, the discriminant (Ik, the bound t{K) we defined above and 
the minimum covering density Q{K). 

Bayer-Fluckiger and Nebe showed by giving a number a, which defines the appro- 
priate inner product (•, that 13 candidates are thin. For computing an upper bound 
of @{K) as given in the table we used their values a G P in all but one case. In the case 
of Q[a;]/(rE'^+x^—3a;—l) their tablecontains amisprint. Instead of a = 1— ISx+lOa;^ 
which does not lie in V we use a = 2 — x instead. For one candidate Bayer-Fluckiger 
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n 


dx 


K 


tiK) 


QiK) 


thin 


2 


5 


0\x]/(x'^ - 5) 

J/ \ / 


> 2.8099 


< 1.2645 


yes 




8 


0\x]/(x'^ - 2) 

J/ V / 


> 2.2214 


< 1.3463 


yes 




12 


0\x]/(x^ - 3) 

J/ v / 

Q[xl/(x2 - 13) 
0\x]/(x'^ - 17) 


> 1.8137 


< 1.2092 


yes 




13 


> 1.7426 


< 1.5708 


yes 




17 


> 1.5238 


< 1.2497 


yes 




21 


X /fx^ - 21) 


> 1.3711 


< 1.2242 


yes 




24 


Ofx /fx^ - 6) 


> 1.2825 


< 1.2583 


yes 


3 


49 


0\x]/(x'^ -\- x'^ — 2x — 1) 


> 3 1093 


< 1.5584 


yes 




81 


Q[x]/(x3 - 3x + 1) 


> 2.4183 


< 2.1225 


yes 




148 


Q[x]/(x3 + x2 - 3x - 1) 


> 1.7891 


< 1.7014 


yes 




169 


Q[x]/(x3 + x2-4x + l) 


< 1.6743 


> 1.8544 


no 


4 


725 


Q[x]/(x* - x-^ - 3x^ + X + 1) 


> 2.9323 


< 2.7045 


yes 




1125 


Q[x]/(x'^ - x^ - 4x2 + 4x + 1) 


> 2.3540 


< 2.2935 


yes 




1600 


Q[x]/(x'^ + 2x^ - 5x2 _ - 1) 


< 1.9740 


> 2.2853 


no 




1957 


Q[x]/(x4 -4x2 -x + 1) 


< 1.7849 


> 1.8939 


no 




2000 


Q[x]/(x'^ - 5x2 _^ 5) 


= 64 


= 64 


weakly 


5 


14641 


Q[x]/(x^ + x'' - 4x^ - 3x^ + 3x + 1) 


> 2.4318 


< 2.2961 


yes 



Table 1 . Classification of totally real thin number fields. 



and Nebe showed that it is weakly thin. In this case t{K) = Q{K) = G4 where 64 
is the least lattice covering density among all 4-dimensional lattices. It is uniquely 
attained by the lattice A|. 

The other three cases were left open by Bayer-Fluckiger and Nebe and by our com- 
putation it turns out that they do not give thin fields. Generally, for a totally real 
number field K, there exists a subspace T of 5" whose dimension equals the degree 
of K, and such that 0(A') = Qt- The corresponding subspace T is given by the basis 

[{uJi^iVj)^^)^^.^.^^, k = 1, . . . ,n, where Or = H + and ai, . . . 

forms a Q-basis of K with G V. Thus if there exist only finitely many inequiva- 
lent T-secondary cones, we can compute G(A') = 0^ by the methods explained in 
Section [6] 

For the field K = Q[x]/ (x'^ + x2 — 4x + 1) we have 7 generic T-secondary cones 
and for all lattices Ax,q our computation proves the bound ©(Aft- ^) > 1.8544. For 
the field K = Q[x]/(x^ + 2x'^ — 5x2 _ g^; — 1) we have 47 generic T-secondary 
cones and for all lattices Kx^a our computation shows 6(Ax,a) > 2.2853. For the 
field K = Q[x]/ (x^ — 4x2 — x + 1) we have 341 generic T-secondary cones and for 
all lattices Aft-,a our computation gives B(Ax,o) > 1.8939. 

The fact that in all these cases the number of inequivalent T-generic secondary 
cones is finite comes as a pleasant surprise. In this situation T is not the space of 
invariant forms of a Bravais group so that we can not a priori rely on the finiteness 
result of Section m It remains an interesting open problem to give a proof of this fact, 
which does not rely on the complete enumeration of T-secondary cones. 
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d 


lattice 


covering density G 


d 


lattice 


covering density G 


1 




1 


13 


1 ^ 

13 


7.762108 


2 


AS 

2 


1.209199 


14 


1 ^ 

14 


8.825210 


3 


AS 


1.463505 


15 


ic 

15 


11.004951 


4 


A! 

4 


1.765529 


16 


A* 

16 


15.310927 


5 


A5 


2.124286 


17 


^17 


12.357468 


6 


■-6 


2.464801 


18 


A* 
^18 


21.840949 


7 


■-7 


2.900024 


19 


AlO 

^19 


21.229200 


8 


■-8 


3.142202 


20 


A^ 
^20 


20.366828 


9 


■-9 


4.268575 


21 


All 
^21 


27.773140 


10 


■-10 


5.154463 


22 


A* 


< 27.8839 


11 


■-11 


5.505591 


23 


A* 


< 15.3218 


12 


■-12 


7.465518 


24 


A24 


7.903536 



Table 2. Least dense known (lattice) coverings up to dimension 24. 



9. Application II: New best known sphere coverings 

In this section we explain how the new theory can be used to construct good and 
even new best known sphere coverings. With the described methods, we were in 
particular able to construct new best known (lattice) sphere coverings in dimensions 
9, . . . , 15. Table 2 gives an overview on the currently best known sphere coverings up 
to dimension 24. This table is an update of the table given in IICS881 Table 2.1]. Note 
in particular, that in comparison to the table there, we have new best known sphere 
coverings in all dimensions d G {6, . . . , 21} \ {16, 18}. Note that the problem has 
been solved for d < 5 only, by the work of Ryshkov and Baranovskii IIRB75II . We 
keep an updated list with additional informations on the involved lattices on our web 
page nSVOSal . 

This section is organized as follows: First we give some background information on 
the lattices of Table 2. Then the next two sections deal with two different constructions 
of subspaces T we used. Both contain specific information on how we obtained the 
new best known lattices. The last section gives some information on a similar approach 
to the closely related packing-covering problem. 



9.1. Some background on best known lattice coverings. The new sphere coverings 
in dimension 6, 7 and 8 were obtained and described in detail in IISV06II and HSVOSbi 
The lattices A^, where r divides d+l, ai^e the Coxeter lattices IICox51ll . A possible 
definition is via the root lattice 

d+l 

Arf = {a;GZ'^+i:^Xi = 0}. 

1=0 

The Coxeter lattice is the lattice generated by A^ and the vector 

/d+l \ 



[l/r) \ - (ei + .... + e(^d+i)/r 



vi=l 
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where denotes the i-th standard basis vector. It is the unique sublattice of con- 
taining J\fi as a sublattice of index r. In particular = A^"*"^. Other well known 
lattices in the series are A^ = E7, Aj = Ey and Ag = Eg. 

Since the symmetric group S^+i acts on the lattice AJJ, it is possible (with the help 
of a computer) to enumerate all its orbits of Delone polytopes and to compute their 
covering densities in fairly large dimensions. A detailed description of the method 
together with results up to dimension 27 can be found in IIDSV07I . 

Baranovskii ||Bar94|| computed (by hand) the Delone decomposition of Ag, finding 
the former best known sphere covering in dimension 9. Anzin IIAnz02l computed 
the Delone decomposition of Af-^ and Ajg establishing the former covering records in 
those dimensions. In a private communication he reported on computing the covering 
densities of A^^ and Afg, which were also best known ones in their dimension. Hence 
suitable Coxeter lattices provide good covering lattices. Nevertheless, by applying 
our theory, we found that all of the five mentioned lattices do not even give a locally 
optimal lattice covering. However, these lattices give good starting points. In fact, we 
obtained the new covering records in dimensions d = 9, 11, 13, 14, 15 by applying our 
theory to a suitable linear subspace T containing a PQF of the corresponding Coxeter 
lattices. More details are given in the next section. For dimensions d = 10, 12 we used 
a lamination technique, which is described thereafter. 

The entries of Table 2 for dimensions d = 22, 23, 24 ai^e "consequences" of the 
existence of the Leech lattice A24. It may not surprise that the Leech lattice itself yields 
the best known lattice covering in dimension 24. Its covering density was computed 
by Conway, Parker and Sloane ( IICS88[ Chapter 23]). In MSVOSbll it is shown that 
the Leech lattice gives at least a local optimum of the covering function Q{Q). Note 
that the root lattice Eg does not have this property. Knowing the comparatively low 
covering density of the Leech lattice. Smith [jSmiSSil was able to estimate the covering 
densities of the dual laminated lattices AS, and AS,. 



9.2. Large subgroups and small linear subspaces. One strategy to find good cover- 
ing lattices is to consider subspaces T containing a PQF, which gives a good or even 
best known sphere covering. In order to keep the number of T-secondary cones low 
(manageable), one preferably chooses a subspace T of low dimension, e.g. less than 
or equal to 4. Another problem is the size of Delone subdivisions to be dealt with. If 
T is contained in a space of invariant forms J^{G) of a preferably large finite subgroup 
G of GL(;(Z), then this problem can be reduced by exploiting the symmetries of the 
Delone subdivisions, respectively those of the forms Q G J^{G) (see Proposition 16.21 ). 

Dimension 9. We optimized over a 3-dimensional subspace T containing a PQF 
of the former record lattice Ag. We computed 210 T-generic T-secondary cones. The 
new covering record is attained by a PQF which has 34 orbits of Delone polytopes. 

Dimension 11. We optimized over a 3-dimensional subspace T containing a PQF of 
the former record lattice Af^. We computed 2444 T-generic T-secondaiy cones. The 
new covering record is attained by a PQF which has 99 orbits of Delone polytopes. 

Dimension 13. We optimized over a 2-dimensional subspace T containing a PQF 
of the former record lattice Ajg. We computed 79 T-generic T-secondary cones. The 
new covering record is attained by a PQF which has 134 orbits of Delone polytopes. 
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Dimension 14. We optimized over a 2-dimensional subspace T containing a PQF 
of the former record lattice /Kf,^. We computed 162 T-generic T-secondary cones. The 
new covering record is attained by a PQF which has 983 orbits of Delone polytopes. 

Dimension 15. We optimized over a 2-dimensional subspace T containing a PQF 
of the former record lattice Afg. We computed 109 T-generic T-secondary cones. The 
new covering record is attained by a PQF which has 203 orbits of Delone polytopes. 

9.3. Best coverings from laminations. Another fruitful strategy is the construction 
of good coverings from lower dimensional ones. Assume that Q G S^q is contained 
in a linear subspace T of S'^. We consider a corresponding lattice L = AU^ with 
A G GL(i(M). Further, we choose a point Ac with c G M*^ and consider all lattices L\ 
in M"'"'"^ generated by 



Choosing Q within T and G M, we obtain a linear- subspace T' of forms Q' with 
dim T' = dim T + 1. If Ac is the center of a Delone polytope P of the lattice L, then 
the lattice vectors in the lattice L\ belonging to the i-th layer 



have at least distance |i|A from (Ac, 0). Hence, P is a Delone polytope embedded in 
L\ when A is at least the circumradius of P in L. 

Dimension 10. We considered the 1 -dimensional linear space T containing Ag. We 
look at the linear subspace T' of S^^ of dimension 2 constructed with the center of a 
Delone polytope having the largest symmetry group in Ag. We computed 6 T'-generic 
T'-secondary cones. The new covering record is attained by a PQF which has 4 orbits 
of Delone polytopes. 

Dimension 12. We considered the 1-dimensional linear space T containing k\^. 
We look at the linear subspace T' of of dimension 2 constructed with the center 
of a Delone polytope. We computed 241 T'-generic T'-secondary cones. The new 
covering record is attained by a PQF which has 1206 orbits of Delone polytopes. Here 
we tried all different centers. The new covering record was produced by using a center 
of a Delone polytope having the third largest symmetry group. 

9.4. New best known packing-coverings. Finally we mention briefly a third appli- 
cation of our new theory. Closely related to the lattice covering problem, is the lattice 
packing-covering problem. It asks to minimize the packing-covering constant 



of a d-dimensional lattice L, where A(L) = ^ min{||?;|| : v G Z \{0}} is the so called 
packing radius. For a detailed description of this problem as well as for its interpre- 
tation as a convex optimization problem, we refer the interested reader to IISV06II . So 




with A > and associated PQF 
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far, optimizing over suitable T-secondary cones with respect to 4-dimensional sub- 
spaces T of 5^, we found a new best known 7-dimensional packing-covering lat- 
tice with 7(L) = 1.499399. . . . The former record holder was the lattice Ey with 
7(Ey) = Y^7/3 = 1.527525. ... It remains to undertake a systematic search for 
further best known, maybe optimal packing-covering lattices. 
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